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Abstract 

Using the conception of the so-cahed quantum relaxation, we build a semiphenomenolog 
ical model for ortho-para conversion of nuclear spin isomers of hydrogen- type molecule. 
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1 Introduction 



Recently an idea was proposed [1] and proved [2-4] that the combined action of intramolecular 
dynamics and environment could account for conversion of nuclear spin modifications of ^^CH^F. 
One may assume the same conversion mechanism for other molecules. This mechanism resides 
in the following. The intramolecular Hamiltonian is nondiagonal with respect to definite total 
nuclear spin J. That implies the nutations between molecular states with various / values. The 
role of the environment is in the braking of nutation's phase. Therefore the environment alone does 
not initiate transitions with nuclear spin value changing, but, destructing quantum coherence, the 
environment brakes the reversibility of nutations. It is remarkable that no cross-section can be 
ascribed to such a conversion way, which makes pointing a classic analogy to this quantum process 
rather problematic. For this reason, we name it quantum relaxation. The known phenomena of 
enantiomers conversion [5], decay of neutral kaons [6], and neutrino oscillations [7] are other 
examples of quantum relaxation. 

Considered as quantum relaxation, the nuclear spin conversion delivers a number of interesting 
problems, especially when being extended beyond the lowest approximations with respect to the 
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spin-mixing intramolecular Hamiltonian. This was shown in the work [8], where quantum relax- 
ation of multilevel system was considered in the context of rather general approach. The price paid 
for the progress, was also rather high. It consisted of specific form of intramolecular Hamiltonian 
which mixed only a pare of states belonging to various spin modifications. The dissipative term 
in the quantum evolution equation, was particular as well, being a modification of BGK collision 
model. 

In the present work an attempt is made to overcome these limitations in the special case 
of hydrogen-type molecule. The success of such an approach is stipulated by specific algebraic 
properties of operators governing the evolution of nuclear spin and molecular frame. The author 
states by no means the quantum relaxation being the main (or even significant) cause for conversion 
of real ortho- and para-hydrogen. We consider if2-type-molecule because of the following three 
reasons: 

1. due to its simplicity if2-nioIccTilc turns to be an ideal testing site for models of the influence 
of the environment and for investigation of sensitivity of the conversion rate to the variety of the 
models; 

2. one can write the most general form for the intramolecular nuclear spin-mixing Hamiltonian 
allowable for H21 

3. it is possible to account and to some extent investigate analitically higher orders of the spin- 
mixing interaction in the conversion process. 

In Sec. 2 and Sec. 3 we introduce the set of operators for, respectively, nuclear spin subsystem 
and molecular frame. In what follows (Sec.4 and below) these operators are used as building 
blocks for construction of spin-mixing Hamiltonian and various dissipative terms. 

2 Operator algebra of nuclear spin subsystem 



The operators acting in the nuclear spin space of H2 molecule are expressed through the pro- 
tons' spin operators /*' and /*" "* and may be considered as a basis of a 16-dimentional associative 
algebra. Let us introduce and describe its elements. This algebra includes naturally the three 

components of the total nuclear spin vector / = /'''* -|- /*■ ^ The total spin magnitude may take 
two values - and 1. One may introduce the operator of the nuclear spin value /: 

/^2+(l'".F>). (2.1) 

with the mentioned eigenvalues. This operator satisfies the evident relation {i ■ 1^ =/(/ + !). We 

will either need raising and lowering operators with respect to /. Not too stretched manipulations 
provide the following vector operators 



with the needed properties, which is evident from the commutators 



[i^^\i] = ^i'^^\ (2.3) 
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Repeated application of any rising and lowering operators' components gives zero: 

jMjM ^ ^ 0. (2.4) 

One has also 
and, respectively, 

fjM^lM^ (2.6) 
The products of raising and lowering operators' components give 



(2.7) 



4 

where a new tensor operator 



4 = |W/f + /<')|f-2i..fi<'>.l<=)) (2.8) 



appears. This operator is at the same time symmetric traceless tensor in 3D space. Supplemented 

by unit, the set {/, /, j''^'*, /*■ \ Pj^} forms the associative algebra of the nuclear spin subsystem 
of E^. 

The following product relations involving P^j will be needed: 











/(-)p.. 


P../(-) - 


■ f A • - n 

■ -'fc -^tj — 





(2.9) 



To finish the set of products which will be used, one should add 

f.f{+) _ j...Ji+) . /(-)/. _ j...J^-) 

r(-) _ f(+) 



/,/]-^ = /r^/,- = 0. (2.10) 



All other operator products are not of interest to us in the present work. 
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3 Operator algebra of molecular frame 



It was shown in [9] that the evolution of molecular rotator could be described within the 

framework of 10-dimentional space spanned by the operator set {J, J , J^^\ J*" "*} which is at the 
same time the Lie algebra of the real simplectic group S'p(4, R) or, speaking more precisely, of its 
simply connected covering metaplectic group Mp(4, M) [10]. Here J is the molecular rotational 

angular momentum; J is the angular momentum operator value, so that {J ■ J) = J(J+1); J 

and J = J ^ are vector operators responsible, respectively, for transitions J J ± 1. 

The structure of the metaplectic Lie algebra is characterized by the following set of commuta- 
tors: 

K±)i 



= -5ij{2J + 1) - 2ieijkJk, 



0; 



(3.1^ 



( + ) 7( + )l 



J<-'1 



(±) 



We need also the rules of J 

{|J,M)}: 



and "* operation in the basis set of molecular rotational states 



J. 



(+) 



{-) 



{-) 



J. 



J,M) 
J,M) 
J,M) 

J,M) 
J,M) 
J,M) 



(J + M + 1)(J + M + 2)|J+1,M + 1), 



{J - M + 1){J - M + 2)\J + 1, M - 1), 



(J + M + 1)(J-M + 1)|J + 1,M), 



(3.2) 



(J- M)(J- M - 1)1 J - 1,M + 1), 



-a/(J + M)(J + M-1)|J - 1,M- 1), 



{J + M){J-M)\J-1,M), 



where J±~^ 



J ]^ i 'lJ'2 J- 



j[ ^ ± i J2 From these relations one can easily prove the 
validity of the following relations which will be widely used below: 



fj X J 



(+)^ 



(J X J 

(y^^y-^) 



=%{.]+ 1) J 

= —tJj , 

J(2J- 1), 
= i(2J- 1)J 



(+) 



J X J) = —tJ J, 



(J^ ^ X J) 

(J^-^ X 



= (J+l)(2J + 3); 
-i(2J + 3)J; 



(3.3) 



/(+) 7(-) _ 72 e I ■ f 7 _ 7 7 . 



(J + i)^5ij - 2(7 + i)tijkJk - Jj- 



J' 



{J ■ r") = {j ■ r ') = 

The last line in ( p.3|) has very clear physical meaning. It states that no scalar operator, which 
acts only in the space of molecular rotation states, can generate transitions with AJ 7^ 0. 
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4 Mixing intramolecular interaction 



For simplicity the model molecular dynamics we are going to consider consists of free molecular 
rotation with a nuclear spin-mixing interaction imposed on. Vibrational and electron degrees of 
freedom are not of interest and the vibrational state is assumed to be symmetric. 

Let us specify the form of the spin-mixing Hamiltonian Hmix- Because of Pauli principle, ortho- 
hydrogen (/ = 1) can have only odd rotational rotational angular momentum values J, whereas 
para-hydrogen (/ = 0) must be in a rotational state with even J. There must be operators in 

Hmix, which generate transitions between ortho- and para-states. I from Sec. 2 are the only 
operators of such a kind. Because of ( p.4|) , the components of these operators may be contained in 
Hmix in the first degree only. The mentioned correlation between / and J should not be violated 

by Hmix- This makes ^ be multiplied by vector operators affecting the molecular rotation by 
changing J by an odd number. We note that any vector operator, which changes J by 3 or more, 

must inevitably contain the scalar products {J^^^ ■ J*"^"*) or (J^ ^ ■ j'^ ■*) and should be equal to 

r(±) . , 1 . , . , ^-(±) 



zero in accordance with (|3.3| ). Hence J are the only operators which / may be multiplied 



by. One can now write the general form of the mixing Hamiltonian: 

Hmix= C.+ (J)(y+^/^^^) + {J^'^ ■i^'^)uJ+{J) 



') + ^ ■i^^^)uj.{j). (4.1) 



The scale factors uj±{J) in ( |4.1|) account for J- dependence of nuclear spin-mixing and can be 



specified by a microscopic model of mixing interaction. The line over the symbols stands for 
complex conjugation. Taken with the rotational Hamiltonian Hq = ujqJ{J + 1), the terms ( |4.1| ) 
add up to the total molecular Hamiltonian H which will be used below. 

A. 

5 Heisenberg operator for nuclear spin value I{t) 



Written in the superoperator form, the kinetic equation for the molecules' density matrix reads 

N 

dtp{t) = cm] ^ ^oim] + E '^nCnimi (5.1) 

ra=l 

where two types of superoperators (Liouvillians) occur: £o[p(i)] = ^i[H,pit)] is the Liouvillian of 
the free Hamiltonian dynamics; £„ are the generators of "irreversible parts" of evolution caused 
by the environment. Various dissipative generators contribute additively and are specified by the 
number n = 1,2, . . . , N; z/„ give the rates of corresponding dissipative processes. In the present 
work we are going to describe three types of dissipative Liouvillians and consider the simplest two 
ones more or less comprehensively. 

In the context of the simplest approach, one assumes that any collision with the environment's 
particles projects the hydrogen molecule onto pure ortho- and para-states; and this is the only 
result of the collision in the considered model. Hence 

c^[p] = ipi + (1 - -i)-p^[i, [p, /]] (5.2) 
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One can easily prove that L\ causes the extinction of ortho-para-states coherence induced by the 
mixing Hamiltonian. 

In contrast to the first model, the second one is formulated in terms of molecular frame 
operators. It postulates that collisions with the environment's particles cause frequent infinitesimal 
rotations of J (but do not affect directly the nuclear spin): 

^[p] = JpJ - + 1)P - \phJ + 1) (5.3) 

Note that the ortho-para-states coherence means at the same time the coherence between states 
with various J (due to the mentioned I — J correlation). Because of the last two terms in (|5.3| ) 
the deorientational collisions destruct the coherence. 

The third dissipative model accounts for collisional transitions changing J. In every collision 
one has J i— J ± 2. So the collisions do not initiate direct transitions between ortho- and para- 
states. The corresponding Liouvillian has the form 

A[p] = [J^^g(J)exp{/3(2J-l)}p,g(J)J^)] (5.4) 
+ [g( J) exp{-/3(2J - 1)} J^^p, J^^\{J)\ + H.c, 

where J^^ = jf^'^ J^^^ are traceless symmetric tensors; /3 = huJo/ksT] q{J) is a function of J. 
This model is a strait analog of that one from the work [9]. The conversion induced by ( |5.5| ) will 
be considered in details elsewhere. 

Arming with Ci and £2, we may return to Eq.( ^.l|) . It is more convenient to deal with 
Heisenberg observable 0{t) rather than with p(t). The corresponding equation of motion for 0(t) 
reads 

dtOit) = C^Oit)], (5.5) 

where the adjoint Liouvillian is introduced. It is determined with respect to the trace scalar 
product Tr{dC[p]) = Tr{C^d]p). Note that Cl = -£o, J^l = jCu and 4 = ^2- 

Of concern to us is the Heisenberg operator I{t) of the nuclear spin value. The detailed analysis 
shows that I(t) (as well as any scalar Heisenberg operator) has the following structure: 

i{t) = A{J, t) + B{J, t)i + c(+)( J, t) ( J^^^ • i^^^) 

+ c^-\j ,t){j^^^ ■ i^'^) + (j^^^ ■ i^'^)c^+\j,t) + (j^^^ • i^^^)&-\j,t) (5.6) 

+D{lt)iJ ■ S) + E{J,t)JiPijJj + F{J,t)jt^Pijfj^^ + j\-^hyj\-^F{lt). 

The coefficients A{J,t), . . . , F{J,t) are to be determined. 

With the formal solution of Eq. (|5.5|) presented as the Taylor expansion 

00 j.n 
n=0 

one may associate the expansions 

CO MU CO j.n 

Ml = E -MJ), • • • , F{J, t) = Y: -Fn{J). (5.8) 

n=0 n=0 
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After rather cumbersome calculations with the use of algebraic relations from Sees. 2 and 3, we 
arrive at the following set of itterative equations: 

A„+i(J) = -4(J + 1)(2J + 3)LW(J + 1) + 4J(2J - 1)L[;\J)- (5.9) 



5„+i(J) = ^J(2J-l)[LW(J)-34-)(J)]-i(J+l)(2J + 3)[4-)(J+l)-3LW (5.10) 



C^uMJ) = -kj + 1)(2^ + 3)^+(J)E„(J) + z(J - 1)(2J - 3)iu4J - 1)F„{J) (5.11) 
o 

+lUJ+{J)[AniJ - 1) - An{J)] - lUJ+{J)Bn{J) + [2iuoqJ - - v^\C\^\j) + «cu+(J)(J + 1)Z^„(J); 



= ^(^ - 1)(2J - 3)^_( J)£;„(J - 1) - %{.] + 1)(2J + 3)cu+(J + 1)F„( J + 1) (5.12) 
o 

-%uj^{J)\An{J - 1) - A„(J)] + %uj-{J)B^{J) + [22a;oJ - ^^i - ^^JC^i^H^) + ^^-{-W - l)Dn{J)] 



D^+,{J) = -2(2J - 1)LW(J) - 2(2J + 3)4-)(J + 1) - i/2/^n(J); (5.13) 

K+i(J) = 4LW(J) - 4Li\j + 1) - 3z/2K(J); (5.14) 

F„+i(J) = [2z(2J - l)c^o - 3z/2]i^n(J) + 4^c^-(J - l)C(+)(^) - 4ia7+( J)C(-)( J - 1), (5.15) 
where the combinations 

lW(J) = -.c.+ (J)CW(J) +zu;+(J)C(+)(J) 

and 

L(-)(J) = -^^_(J)Ci-)(J) +za;_(J)C(-)(J) 

are introduced. In the case of I{t) we must use the initial conditions when all coefficients 
Ao{J), . . ., Fo{J) are zero except Bo{J) = 1. In principle, Eqs. (|5.9| - |5.15|) let one evaluate I(t) 
provided the Taylor expansions converge. 
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6 Conclusion 



Fundamentally, the Heisenberg operator i(t) let one evaluate any nuclear-spin-dependent mag- 
nitude. For example, the expression 

(/(t)) = Tr/(t)p (6.1) 

gives the conversion process in an initially prepared non-equilibrium state, where p is the initial 
density matrix. One can easily show that (due to the I-J correlations) any relevant p must obey 
the equality 

(2/ + (-1)^) p = p = p{2i+ i-iy) (6.2) 

The stationary density matrix annulates both the left- and right-hand sides of ( p.l| ). It has a 
form akin to (|5l6| ) with t-independent coefficients. The simplest acceptable (but not stationary) 
isotropic density matrix is 

p = PoiJ)+piiJ)i (6.3) 

where in accordance with ( |6.2| ) po{J) =0 for odd J and po{J) = —pi{J) for even J. 

The proposed model allows rather extended analitical evaluations. Numerical calculations are 
still inevitable. For this one should specify the factors ujj-{J) in ( |4.1| ). This is an independent 
problem. There can be faced a situation of bad convergency of the Taylor expansions ( ^.9|) for 
some uj±{J). The analysis of these and related problems will be done elsewhere. As a result we 
will be able to estimate the value of the model completely. 
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